Experimental verification of a modified fluctuation-dissipation relation for a 
micron-sized particle in a non-equilibrium steady state 
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A modified fluctuation-dissipation-theorem (MFDT) for a non- equilibrium steady state (NESS) 
is experimentally checked by studying the position fluctuations of a colloidal particle whose motion 
is confined in a toroidal optical trap. The NESS is generated by means of a rotating laser beam 
which exerts on the particle a sinusoidal conservative force plus a constant non-conservative one. 
The MFDT is shown to be perfectly verified by the experimental data. It can be interpreted as an 
equilibrium-like fluctuation-dissipation relation in the Lagrangian frame of the mean local velocity 
of the particle. 



X 



The validity of the fluctuation-dissipation theorem 
(FDT) in systems out of thermal equilibrium has been 
the subject of intensive study during the last years. We 
recall that for a system in equilibrium with a thermal 
bath at temperature T the FDT establishes a simple re- 
lation between the 2-time correlation function C(t — s) 
of a given observable and the linear response function 
R(t — s) of this observable to a weak external perturba- 
tion 



MFDT reads for t > s 



d s C(t- s) = k B TR(t- s). 



(1) 



However, Eq. ([]} is not necessarily fulfilled out of equi- 
librium and violations are observed in a variety of sys- 
tems such as glassy materials [l|, H, H, fl HI, granu- 
lar matter and biophysical systems Q. This mo- 
tivated a theoretical work devoted to a search of a gen- 
eral framework describ ing FD relations, see the review 
H or @, [H GH, GJ, H G3 for recent attempts in 
simple stochastic systems. In the same spirit, a modi- 
fied fluctuation-dissipation theorem (MFDT) has been 
recently formulated for a non-equilibrium steady dy- 
namics governed b y t he Langevin equation with non- 
conservative forces [la]. In particular, this MFDT holds 
for the overdamped motion of a particle on a circle, with 
angular position 9, in the presence of a periodic poten- 
tial H(9) = H(6 + 2ir) and a constant non-conservative 
force F 



d=-d H(d) + F + (, 



(2) 



where £ is a white noise term of mean (£ t ) = and 
covariance (CtCs) = 2DS(t — s), with D the (bare) diffu- 
sivity. This is a system that may exhibit an increase in 
the effective diffusivity |T^, |T3] - Here, we shall study the 
dynamical non-equilibrium steady state (NESS) reached 
for observables that depend only on the particle posi- 
tion on the circle so are periodic functions of the angle 
9. Such a state corresponds to a constant non- vanishing 
probability current j along the circle and a periodic in- 
variant probability density function po(9) that allow us 
to define a mean local velocity vq(9) = j/po(9). This is 
the average velocity of the particle at 9. For a stochas- 
tic system in NESS evolving according to Eq. ©, the 



d s C(t - s) - b(t - s) = k B TR{t - s), 



(3) 



where the 2-time correlation of a given observable 0(9) 
is defined by 



C(t-s) = (O(9 t )O(9 s )) , 



(4) 



and the linear response function to a 5-perturbation of 
the conjugated variable h t is given by the functional 
derivative 



[0(9 t )) f 



(5) 



h=0 



In Eq. ([5]), (•••)/» denotes the average in the perturbed 
time-dependent state obtained from the NESS by re- 
placing H{9) in Eq. © by H{9) - h t O{9). It reduces 
for h = to the NESS average (...) . In Eq. ©, the 
correlation b(t — s) is given by 



b(t-s) = (0(9 t )v (9 s )deO(9 s )}o. 



(6) 



This new term takes into account the extent of the vi- 
olation of the usual fluctuation-dissipation relation (|TJ) 
due to the probability current and it plays the role of a 
corrective term to C(t — s) in the MFDT, Eq. ([3|), which 
can be rewritten in the integral form: 



C(0) - C(t) - B(t) = k B T X (t), 



(7) 



where B(t) = J* * b(t — s)ds and x(i) = f Q R(t — s)ds is 
the integrated response function. 

In this letter, we present an experimental test of 
Eq. in the linear response regime around a NESS 
attained by a micron-sized particle in a toroidal opti- 
cal trap similar to the one used in [ijj. We first show 
that the dynamics of the particle is well described by the 
Langevin equation © on a circle. Secondly, by measur- 
ing vo, B(t), C(t) and x(i), we verify Eq. (J7|) for the 
observable 0(9) = sm.9. The result can be interpreted 
as an equilibrium-like fluctuation-dissipation relation in 
the Lagrangian frame of the mean local velocity ^o(^) 
[l5[- We also check that po{9) is frame invariant. 
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Figure 1: (a) Invariant density of the angular position of the 
particle in NESS for a modulation of 7% (black solid line) 
and 7.35% (red dashed line) around the laser mean power. 
Inset: Mean local velocity of the particle in the former case, 
(b) Corresponding potential profiles. The arrow indicates 
the direction of the non-conservative force /. 



The experiment is performed using a spherical silica 
particle of radius r =1 pica in ultrapure water at room 
temperature T = 20.0 ± 0.5°C at which the dynamic 
viscosity of water is r\ = (1.002 =p 0.010) x 10~ 3 Pa s. 
The particle is kept by an optical tweezers in a toroidal 
optical trap. This kind of trap consists on a Nd:YAG 
diode pumped solid state laser beam (A = 1064 nm) 
which is focused by a microscope objective (63 x, NA 
= 1.4) and scans (by means of two acousto-optic deflec- 
tors) a circle of radius a = 4.12 pm in the horizontal 
plane at a rotation frequency of 200 Hz. The toroidal 
trap is created 10 pm above the inner bottom surface of 
the cell where hydrodynamic boundary-coupling effects 
on the particle are negligible. At a rotation frequency 
of 200 Hz, the laser beam is not able to hold the parti- 
cle but drags it regularly a small distance on the circle 
when passing through it [l8[. The diffusive motion of 
the particle along the radial and vertical directions dur- 
ing the absence of the beam is less than 40 nm, thus 
the angular position of the particle 9 (measured modulo 
2-71") is the only relevant degree of freedom of the dynam- 
ics. The laser power is sinusoidally modulated around 
30 mW with an amplitude of 7% of the mean power, 
synchronously with the deflection of the beam at 200 
Hz creating a static sinusoidal intensity profile along the 
circle. This trapping situation acts as a constant non- 
conservative force / associated to the mean kick which 
drives the particle across a sinusoidal potential U (9) due 
to the periodic intensity profile. The particle barycenter 
(xt, yt) is measured by image analysis with an accuracy 



of the order of 1 nm at a sampling rate of 150 Hz and 
exposure time of 1/300 s. This measure allows us to 
determine the angular position of the particle 9t with 
respect to the trap center. For more details about the 
experimental apparatus see Ref. We determine the 
value of / and the profile oiU{9) by means of the method 
described in [2(|. This method exploits the probability 
current j and the invariant density po(9) in NESS to 
reconstruct the actual energy landscape of the particle. 
We recorded 200 time series {9 t } of duration 66.67 s 
with different initial conditions {9q} sampled every 5 
minutes in order to measure j and p$ (9) . The probabil- 
ity current is related to the global mean velocity of the 
particle by the expression j = (9)o/(2ir). The value of 
(9)o is calculated from the slope of the linear fit of the 
mean angular position of the particle (not taken modulo 
27r) as a function of time leading to j = 3.76 x 10~ 2 s _1 
in the direction of the laser beam rotation. The invari- 
ant density, shown as a solid black line in Fig. QJa), is 
computed from the histogram of each time series {9 t } 
averaged over the 200 different initial conditions. In 
Fig. [3Ja) we & lso show the mean local velocity vq(9) = 
j/po(9) of the particle. From these quantities we obtain 
/ = 3r F aj f£* p (9')~ 1 d9' = 6.60 x 10" 14 N and U{9) = 

~k B T\ogp (9) + J°(f - 6irr 1 rajp (9')- 1 )ad9' = Asin(9 
with amplitude A = 68.&ksT. The experimental poten- 
tial profile is shown in Fig.QJb) (black solid line). Hence, 
the time evolution of 9 is claimed to follow the Langevin 
dynamics of Eq. © [H with F = f/(6nr)ra) = 0.85 
rad s" 1 , H{9) = U{9)/{&Tn 1 ra 2 ) = H sin9, H = 
A/(6m]ra 2 ) = 0.87 rad s" 1 , and D = k B T/{Qm ir a 2 ) = 
1.26 x 10 -2 rad 2 s" 1 . Note that in the corresponding 
equilibrium situation (/ = 0) the probability maximum 
would be located at the minimum of U(9) (9 = 3ir/2). 
However, in NESS the non-conservative force / > 
shifts the maximum of po(9) in the positive direction, 
as shown in Fig. Ufa). The choice of the parameters 
has been done to enhance the stochastic nature of the 
dynamics, i.e. we take (F — H )/H n w —2% which is 
close to the maximum increase of the effective diffusivity 

following Rcfs. mm- 

Additionally, 500 times series of duration 100 s were 
specially devoted for the determination of x(f). In this 
case, during each interval of 100 s we apply from time 
to to t + At with < i < 66.67 s and At = 33.33 s 
a step perturbation changing the value of A to A + 5 A. 
This is accomplished by suddenly switching the laser 
power modulation from 7% to 7.35% of the mean power 
(30 mW) . The experimental value of the amplitude per- 
turbation (5 A = 0.05 A) is determined from independent 
NESS measurements of po(9) and U (9) for a power mod- 
ulation of 7.35% (shown in Figs. H^a) and HJb) respec- 
tively as red dashed lines) as described previously. By 
keeping constant the mean power during the switch we 
ensure that the value of / remains also constant, com- 
pare to a different time-dependent protocol explored re- 
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Figure 2: (a) Comparison between the different terms needed 
to verify Eq. ||7J, as functions of the time lag t. (b) Expanded 
view of the comparison between the curves C(0) — C(t) — B(t) 
and ksTxit) shown in Fig[2ja). The thin red dashed lines 
represent the error bars of the measurement of the integrated 
response. 




Figure 3: Example of trajectories used to compute the in- 
tegrated response. We show a perturbed trajectory (thick 
dashed red line) and four out of a total of 200 of the corre- 
sponding unperturbed ones (see text). 



cently in [2l| . In this way, we extract 500 perturbed 
trajectories {9 t }sA of duration At = 33.33 s. We ensure 
that after switching off the perturbation the system ac- 
tually has attained a NESS before the beginning of the 
next step perturbation. 

With the purpose of determining correctly the differ- 
ent terms involved in Eq. ([7]), the observable 0(9) must 
be chosen consistently on both sides of such relation. 
The change of the potential U(9) -> U(9) + 6Asm9 due 
to the application of the step perturbation implies that 
0{9) = sin (9 is the observable that must be studied with 
— SA as its conjugate variable. Hence, we compute the 
correlation function C{t), the corrective term B{t) (with 
dgO{9) = cos 9 and the experimental curve vq(9) shown 



in Fig. Ufa)) and the integrated response x(t) for this 
observable, as functions of the time lag t. 

The determination of C(t) and B(t) is straightforward 
according to Eqs. (U]) and ©■ The stationarity of the 
system allows us to perform an average over the time 
origin in addition to the ensemble average (...)o over the 
200 different time series devoted to this purpose, which 
increases enormously the statistics. The dependence of 
C(0) — C{t) and B(t) on t is shown in Fig. [UJa) in green 
dotted-dashed and blue dashed lines, respectively. 

On the other hand, the integrated response x(t) is 
given by 



x(t) 



(0(9t))sA-(0(9t))o 
-SA 



(8) 



In Eq. ([8]) the value t — corresponds to instant to 
when the perturbation of the potential amplitude 5 A is 
switched on. To decrease the statistical errors in com- 
parison of the terms in Eq. ([H]), for a given perturbed 
trajectory 9tsA we look for as many unperturbed ones 
9t as possible among the 200 time series {9t} starting at 
a time t* such that 0(9t*) = 0{9osa)- Then we redefine 
t* as t = in Eq. (JH), as shown in Fig. [3l The unper- 
turbed trajectories found in this way allow us to define a 
subensemble over which the average (O(9 t )) in Eq. © 
is computed at a given t. The average {0(9 t ))sA is sim- 
ply computed over the 500 perturbed time series. In Fig. 
HJa) we show in thick dashed red line the dependence of 
the integrated response on t. 

The comparison between the different terms needed 
to verify Eq. ((3|) is shown in Fig. [2ja), for the time lag 
interval < t < 3.5 s. As expected, the usual FD re- 
lation (JT|) is strongly violated in this NESS because of 
the broken detailed balance, with the correlation term 
C(0) — C(t) being one order of magnitude larger than 
the response term ksTxit)- However, with the correc- 
tive term B(t) associated to the probability current sub- 
tracted, C(0) —C{t) — B(t) shown in solid black line in 
Fig.[UJa), becomes equal to ksTxit)- For clarity, in Fig. 
[3Jb) we show an expanded view of the of the curves 
C(0) - C{t) - B(t) and k B T X (t). We observe that, 
within the experimental error bars, the agreement be- 
tween both terms is quite good, verifying the integrated 
form of the modified FD relation ([7]). The error bars 
of the integrated response curve at each time lag t arc 
obtained from the standard deviation of the subensem- 
ble of unperturbed trajectories found for each perturbed 
trajectory, like the ones shown in thin solid lines in 
Fig. El We checked that the perturbation SA = 0.05A 
is small enough to remain within the linear response 
regime. This is quantitatively seen in Fig. 2(b) where 
the response measured at SA = 0.0 7 A is represented by 
circles. We see that %(£) is independent of SA within 
experimental errors showing that we are in the linear 
response regime. The MFDT is checked only for the 
first 3.5 s because after this time the evaluation of x(i) 
is affected by large errors. Indeed during the measure- 
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Figure 4: Invariant density of the angular position measured 
in the Eulerian frame (continuous line) (same as in Fig.fjja)) 
and in the Lagrangian frame (*). Inset : example of a tra- 
jectory measured respectively in the Eulerian (a) and La- 
grangian (b) frames. 

ment of x(*)> f° r finite 8 A the system is approaching a 
new NESS which depends non-linearly on 6 A (see the 
strong non-linear dependence of po on 6 A in Fig.[TJa)). 
The perturbed trajectories diverge with respect to the 
unperturbed ones (Fig. 3) in such a way that after 3.5 
s the non-linear effects and the statistical error bars of 
\(t) due to finite sampling become large. 

As shown in the validity of Eq. Q for the fluc- 
tuations of the angular position of the silica particle in 
NESS gains a simpler interpretation in the Lagrangian 
frame of the mean local velocity vq(9) along the circle. 

Indeed, using the observables that are time indepen- 



dent in the Lagrangian frame, the MFDT may be rewrit- 
ten in the form 

d s C L (t,s) =k B TR L (t,s), (9) 

where Cl and Rl are the correlation and the response 
measured in the Lagrangian frame. Eq. ^ is close to 
that of the equilibrium FDT (Eq. (fT|)) except for the 
lack of the time translation invariance of the functions 
involved. One of the new predictions of the Lagrangian 
analysis of the system is that, although the trajectories 
in the Eulerian and the Lagrangian frame are quite dif- 
ferent, their average density po is the same in the two 
frames. This property is clearly illustrated by the exper- 
imental data in Fig. 01 where we compare the densities 
measured in the two frames. The insets of Fig. 0] point 
out to the difference between a trajectory measured in 
the Eulerian frame and the same trajectory measured in 
the Lagrangian frame. 

We have verified experimentally a modified 
fluctuation-dissipation relation describing the dy- 
namics of a system with one degree of freedom in NESS, 
namely a Brownian particle moving in a toroidal optical 
trap. We point out that the experimental results 
reported here represent an alternative approach to 
non-equilibrium fluctuation-dissipation relations to that 
of Ref. [l^l which dealt with the velocity fluctuations 
relative to the mean local velocity. The approach 
followed in our work relies on an observable depending 
on the particle position. It quantifies the extent of the 
violation of the usual FDT by means of the term B(t), 
admitting a transparent Lagrangian interpretation of 
the resulting MFDT. 
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